PROCEEDINGS 

OF SCIENCE 



Nucleon form factors and moments of parton 
distributions in twisted mass lattice QCD 



o 

(N 
o 

Q 



(N 
O 

q 
o 

(N 



C. Alexandrou {a < b \ M. Constantinou* {a \ C. Kallidonis W, T. Korzec ^ 

^ Department of Physics, University of Cyprus, P.O. Box 20537, 1678 Nicosia, Cyprus 

W Computation-based Science and Technology Research Center, Cyprus Institute, 20 Kavafi 

Str., Nicosia 2121, Cyprus 

( e ) Institut fur Physik Humboldt Universitat zu Berlin, Newtonstrasse 15, 12489 Berlin, Germany 



E-mail: alexand@ucy.ac.cy constant inou . martha@ucy . ac . cy 



kallidonis .christos@ucy.ac.cy korzec@physik . hu-berlin . de 



M. Brinet, J. Carbonell, P. A. Harraud, M. Papinutto 

Laboratoire de Physique Subatomique et Cosmologie, UJF/CNRS/IN2P3, 53 avenue des 
Martyrs, 38026 Grenoble, France 



£ , -mfli7:nariane@lpsc . in2p3 . f r Jaume . CarbonellS Ipse . in2p3 . f r 



harraud@lpsc . in2p3 . f r Ylauro . Papinutto@ Ipse . in2p3 . f r 



P. Guichon 

CEA-Sa clay IRFU/SPhN, 91191 Gif-sur -Yvette, France 
E-mail: |pierre . guichongcea . f r 



K. Jansen 

NIC, DESY, Platanenallee 6, D-15738 Zeuthen, Germany 



E-mail: Karl.Jansen@desy.de 



We present results on the electroweak form factors and on the lower moments of parton distri- 
butions of the nucleon, within lattice QCD using two dynamical flavors of degenerate twisted 
mass fermions. Results are obtained on lattices with three different values of the lattice spacings, 
namely a=0.089 fm, a=0.070 fm and a=0.056 fm, allowing the investigation of cut-off effects. The 
volume dependence is examined by comparing results on two lattices of spatial length L=2.1 fm 
and L=2.8 fm. The simulations span pion masses in the range of 260-470 MeV. Our results are 
renormalized non-perturbatively and the values are given in the MS-scheme at a scale pi=2 GeV. 
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1. Introduction 

Lattice QCD calculations of observables related to the structure of baryons are now being 
carried out using simulations of the theory with parameters close to the physical regime both 
in terms of pion mass and with respect to the continuum limit [||, g[ |[ |], [50. In particular, a 
number of major collaborations are investigating the form factors and the first moments of par- 
ton distributions of the nucleon that encode important information on its structure. In this work 

ill U - III ^Mn} 

we consider the nucleon matrix elements of the operators Vy = YY iD ...iD \j/ and 

v£ = \j/y^ l iD . . .iD 75 \j/ limiting ourselves to n = 1 and n = 2. The case n = 1 reduces 
to the nucleon form factors. The nucleon matrix element of the electromagnetic current is written 
in the form 

(N(p>,s>) \/\N(p,s)) = ( ^ V 2 %(p y) 



Y^{q 2 )+'^—^F 2 ( 
2niN 



u N (p,s) , 
(1.1) 

where q 2 = (p' — p) 2 , is the nucleon mass and £/v(p) its energy. F\ (0) = 1 for the proton and 
zero for the neutron since we have a conserved current. F 2 (0) measures the anomalous magnetic 

moment. They are connected to the electric, Ge, and magnetic, Gm, Sachs form factors by the 

relations: Gpiq 2 ) = F\{q 2 ) + ^ N ^ F 2 (q 2 ) anc * ^m(<7 2 ) = F\(q 2 ) +F2(q 2 ). Similarly the matrix 

element for the axial-vector current is given by 

,2 \ V2 



(N(p',s')\Al\N(p,s))=i[ N u N (p',s') 



G A (q 2 )Y,Y5 + P^G p (q 2 ) 



-u N (p,s). 



(1.2) 

We also study matrix elements of operators with a single derivative (n = 2). The matrix elements of 
these operators are parameterized in terms of the generalized form factors A2o(G 2 ), fi2o(G 2 ) 5 C2()(G 2 ) 

and A 20 (G 2 ), B 20 {Q 2 ) (Q 2 = -q 2 ). The results are non-perturbatively renormalized [[J. In this 
work we use simulations with Np =2 twisted mass fermions (TMF). They provide an attractive 
formulation of lattice QCD that allows for automatic €?(a) improvement, infrared regularization 
of small eigenvalues and fast dynamical simulations [[J. For the calculation of the nucleon ob- 
servables discussed in this work the automatic &{a) improvement is particularly relevant since it 
is achieved by tuning only one parameter in the action, requiring no further improvements on the 
operator level. Details on the lattice evaluation can be found in Refs. [g]. 

2. Nucleon form factors 

Our lattice results on the nucleon axial charge, gA are shown in Fig. |I[ We observe that results 
at our three different lattice spacings are within error bars and that results at the two different 
volumes are also consistent. In the same figure we compare to the results obtained using Np = 
2 + 1 domain wall fermions (DWF) by the RBC-UKQCD collaborations [Q] and using a hybrid 
action with 2+1 flavors of staggered sea and domain wall valence fermions by LHPC [Jsj] . There 
is agreement among lattice results using different lattice actions even before taking the continuum 
and infinite volume limit indicating that lattice artifacts are small. The pion mass dependence for 
the nucleon axial charge has been studied within HB^PT in the SSE formulation [St]. Fitting the 
volume corrected and continuum extrapolated TMF results shown in Fig. [j] we find £4 =1.12(8), 
while fitting the raw lattice results we obtained g&= 1.08(8) confirming that cut-off and volume 
corrections have indeed a small effect. 
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Figure 1: Left: Results on gA as a function of the pion mass squared using TMF, DWF and domain 
wall valence on a staggered sea (hybrid). Right: gA obtained by taking the continuum limit of the volume 
corrected TMF data. The shaded area is the best chiral fit to the data shown on the graph. 



1.0 
0.8 

CM 

3 0.6 

Q. 

o" 0.4 
0.2 
0.0 



V 



■ N,=2 TMF a=0.089fm L=2.8fm 
X N,=2+l DWF a=0.084 L=2.7fm 

° N p =2+1 Hybrid a=0.124 L=2.5fm 
* N =2 Clover a=0.07 L=3.4fm 



I 



' experiment 



1.2 

1.0 

g 0.8 

o" 0.6 
0.4 
0.2 



N =2 TMF a=0.089fm L=2.8fm m =298 MeV 



X N =2+1 DWF a=0.084fm L=2.7fm m =330 MeV 



Q 2 (GeV 2 ) 



1.0 1.5 

Q 2 (GeV 2 ) 



\ 



■ N f =2 TMF a=0.08Bfm L=2.8fm 
X N f =2+1 DWF a=0.084 L=2.7fm 

O N F =2+1 Hybrid a=0.124 L=2.5fm 
* N =2 Qover a=0.07 L=3.4fm 



^ 15 

s 

tD^ 10 



■ N F =2 TMF a=0.0B9fm L=2.Bfm 111^=298 MeV 

A N =2+1 Hybrid a=0.124 L=3.5fm m =353 MeV 



0.6 0.8 1.0 

Q 2 (GeV 2 ) 




1.5 

Q 2 (GeV 2 ) 



Figure 2: Ge and Gm at m K ~300 MeV (left) and Ga and G p (right) for various lattice actions. The filled 
(green) circles show experimental data on Ge and Gm- The solid line on the right panel is a dipole fit to 
experimental data for Ga{Q 2 ) combined with pion pole dominance to get the solid curve shown for G P (Q 2 ). 



In Fig. U we show the isovector electric and magnetic form factors at m n ~ 300 MeV as a 
function of the momentum transfer squared. We compare our results with results using other 6(a) 
improved actions at the same pion mass [|8 , |l0 f 11]. The results are in agreement in the case of Ge, 
while in the case of Gm clover results from Ref. [ 1C] are lower. In the same figure we also show 
TMF results on the two axial- vector form factors and we compare with DWF at m n = 330 MeV [Q] 
and hybrid results by LHPC at m % = 356 MeV ^ on a lattice with L = 3.5 fm. The results are in 
agreement in the case 

°f Ga(Q 2 ), while in the case of G P {Q 2 ) there are discrepancies at low Q 2 
values, which may indicate that volume effects are not negligible on form factors such as G P (Q 2 ) 
which are strongly affected by the pion-pole. 
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Figure 3: Our lattice data on (x)„_ f / (A20) and (x)a„-ac/ (A 20) as a function of the pion mass squared. Results 
using other lattice actions are also plotted. The physical point is shown by the asterisk. 

3. Nucleon moments 

In Fig. H] we show our results on the spin-independent and helicity moments, and results using 
other lattice actions [|L 12, [13|]. Although to compare lattice data using different discretization 
schemes one would have to first extrapolate to the continuum limit, we find a good agreement 
among lattice results, due to small cut-off effects for lattice spacings of about 0.1 fm. Lattice 
values for (x) H _j = A2o(<2 2 = 0) although compatible, are higher from the phenomenological value 
{x u -d) ~ 0.16. A similar conclusion holds for the helicity moment. 
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